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ABSTRACT 

o : 

We present a theoretical analysis of Rossi X-ray Timing Explorer data of Z source 
GX 340+0 obtained by Jonker et al. In the frameworks of the recently formulated the 
transition layer model the 5— angle is an angle between the neutron star (NS) magne- 

^ ■ tospheric axis and the disk (presumably NS rotational) axis. We determine the angle 

Q\ • 5 = 6°. 3 ± 0.°3 which is a combination of the simultaneously observed kHz QPO and 

HBO frequencies. While these three frequencies change by a factor of three or more 

^—i . their 5— combination stays almost constant. GX 340+0 is the fourth source (in addi- 

tion to 4U 0614+09, Sco X-l and 4U 1702-42) for which 5 has been determined. With 
at most one (constrained) parameter we make a complete classification of six observed 
power spectral features, including the two kHz frequencies, the first and second harmon- 
ics of the HBO frequency, low-frequency noise component and break frequencies. We 
demonstrate that a new component discovered by Jonker et al. in the GX 340+0 power 

-»». ■ spectrum is related to the viscous frequency branch which has been, in fact reported 

earlier in 4U 1728-34 by Ford and van der Klis (1998). Finally, we re-classify several 
q \ previously misidentified features in the power spectrum. 

a: 

Subject headings: accretion, accretion disks — diffusion — stars:individual (GX 340+0, 



4U 0614+09, 4U 1728-34, Sco X-l, 4U 1702-42)— stars:neutron— X-ray:star— waves 



Introduction 



This Letter contains the classification of power spectral features in GX 340+0 and a determi- 
nation of the 5— invariant using observational data from Jonker et al. (2000). Psaltis, Belloni & 
van der Klis (1999) discussed many similarities and correlations between the QPO frequencies in 
atoll sources, Z-sources and black hole candidates, which may be related to the similar phenom- 
ena and physical processes occurring in all of these systems. In addition, Wijnands & van der Klis 
(1999) found that the break frequency of the broken power law which describes the power spectrum 
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correlates well with the frequencies of the peaked noise components observed in atoll sources and 
black hole candidates. Titarchuk & Osherovich (1999; hereafter T099) attempted to classify QPO 
frequencies using the basic properties of accretion flow around compact objects. The first property 
is an adjustment of the Keplerian flow in the disk to the sub-Keplerian inner boundary condition 
(Titarchuk, Lapidus & Muslimov 1998, hereafter TLM98), which is a common phenomenon for 
neutron star (NS) and black hole candidate systems. In this respect all correlations between the 
break frequencies Vb and low QPO frequencies vy should be similar in both types of system be- 
cause they are related to the diffusion and oscillation time scales in the transition region between 
the Keplerian flow and the rotational inner boundary (either the NS surface or the inner edge of 
the accretion flow in BH systems). Later the predicted correlation was confirmed by Ford & van 
der Klis (1998) who analyzed RXTE data from 4U 1728-32. Furthermore, in T099 the authors 
demonstrated that and vy are related through the power law oc Uy & which fits to the data in 
most BH and NS systems (Wijnands Sz van der Klis, 1999). 

The second property is an effect of the sub-Keplerian rotating magnetosphere located above 
the disk. Osherovich & Titarchuk (1999a, hereafter OT99a) suggested that the radial oscillations of 
the fluid element bounced from the disk shock region (presumably at the adjustment radius) would 
be seen as two independent oscillations in the radial and the vertical directions due to the presence 
of Coriolis force in the magnetospheric rotational frame of reference. The observed evolution of two 
kHz peaks allows one to infer the profile of the differentially rotated magnetosphere and then to 
verify it using the observed correlation of the horizontal branch quasi-periodic oscillations (HBO) 
frequencies with the low kHz peak. 

In a black hole system the corona above the disk should also rotate with less than Keplerian 
velocity. When oscillating fluid elements of the disk are thrown to the corona, these oscillators 
would be under the influence of the Coriolis force which leads to the same rotational splitting effect 
(OT99a) seen in the NS systems. The main QPO Keplerian frequency would be split into two 
frequencies, the hybrid and low branch related to the radial and vertical oscillations respectively in 
the rotational frame of references. 

The recent discovery of a 450 Hz QPO (Strohmayer 2001) in addition to the previously reported 
18 and 300 Hz (Remmilard 1999) in GRO J 1655-40 provides additional support for the rotational 
splitting effect. 

There are three other QPO models in the literature. The first one is the beat frequency model 
suggested by Alpar & Shaham (1985), and recently modified as a sonic-point beat frequency model 
by Miller, Lamb & Psaltis (1998). The second is the relativistic precession model proposed by 
Stella & Vietri (1998, 1999). The third is the inner accretion disk model promoted by Psaltis & 
Norman (2000). All observational arguments for these models and their problems are discussed in 
detail by Wu (2001) where these models are also compared with the transition layer model (TLM) 
(TLM98, T099-OT99a). We do not go into the details of this comparison in this Letter but we do 
give more arguments for TLM in terms of the new QPO data for GX 340+0 (Jonker et al. 2000, 
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here after JOO). 

Precise simultaneous measurements of the frequencies of the two kHz QPOs and HBO harmon- 
ics in a wide frequency range for GX 340+0 (JOO) give us an opportunity to derive the 5— invariant. 
In Osherovich & Titarchuk (1999, hereafter OT99b), the authors demonstrated for the source 4U 
1702-42 that the inferred angle 5, (see Eq 2, 5 in OT99b) 



S = arcsin 



(^-4r 1/2 (*w^)l (i) 



depends only on observed frequencies: the low and high kHz QPO peaks z/k, and the HBO (the 
low branch) peak vl. It stays the same (3.9° ± 0.2°) over significant range of z^k (650-900 Hz). 
Using a series of observations where v^, uk, vl are detected simultaneously [van der Klis et al. 
1997; Markwardt, Strohmayer & Swank 1999; van Straaten et al. 2000; J00)] we check that the 
angle 5 is a true invariant for four specific sources: Sco X-l, 4U 1702-42, 4U 0614+09 and GX 
340+0 (see §2). We present the magnetospheric rotational profile for GX 340+0 inferred from JOO 
data in §3. Using this profile and the 5 value we will reproduce the theoretical dependence of the 
low branch on the low kHz frequency in §4 and compare it with the JOO observations. 

According to the TLM, also known as "two-oscillator model" all frequencies (namely v^, vl, 
and i/y) have specific dependences on ^k- 

The comprehensive summary of the TML model and its successes to date is present in Titarchuk 
k Osherovich (2000a) 

In §5 we offer the full classification of the QPO frequencies containing six branches which fit 
to the data with at most one (constrained) parameter. Discussion and summary follow in the last 
section. 



2. 5— Invariant and Verification of Transition Layer Model 

We have used the frequencies measured for sources GX 340+0, 4U 0614+09, Sco X-l and 4U 
1702-42 where the low and high kHz QPO peaks vk and and HBO frequencies vl are measured 
simultaneously. The resulting values of S calculated from Eq. (1) are shown in Figure 1. Indeed, 
for each of these sources, the 5— values show little variation with j/k, Vh, V L- They are 3°. 9 ± 0°.2, 
5°. 5 + 0°. 5, 15°. 3 + 0°. 5, for 4U 1702-42, Sco X-l and 4U 0614+09 respectively. The angle 5 obtained 
for the source GX 340+0 is 6°. 3 ± 0°.3 and is similar to the values obtained for Sco X-l and 4U 
1702-42. The existence of the 5— invariant predicted by the TLM [OT99b, Titarchuk & Osherovich 
(2000b here after TO00)] is a challenge for any other QPO model. It is important to note that all 
frequencies included in this S— relation are observed frequencies and thus the relation Eq. (1) is a 
model independent invariant. It is worth pointing out once again that the delta— angle varies from 
source to source. 

In fact, Wu (2001) confirms the existence of the delta-invariant for three more sources - Cyg 
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X-2, GX 17+2 and GX 5-1, using the RXTE observational data obtained by Psaltis et al. (1999), 
Wijnands et al. ( 1998a, b, 1997). For all these sources the 5— angle is about 6° which is very close 
to that obtained for Sco X-l and GX 340+0. The accuracy of 5— determination is of order 5%. 



3. Inferred Rotational Frequency Profile of the NS Magnetosphere in GX 340+0 

From the observed kHz frequencies (J00, Table 3), and z^k, the profile of v mag = Vl{vk)/2i{ 
has been calculated according to 

^ = [4 + 4(^/2tt) 2 ] 1/2 , (2) 

and modeled using the theoretically inferred magnetic multipole structure of a differentially rotating 
magnetosphere (OT99a) 

Vmag = C + C\V^ + C 2 V^[ 3 + C^ K (3) 

where C 2 = 2(C 1 C 3 ) 1 / 2 . The constants C = v mag = 330 Hz, C x = -9.72 x 10~ 2 Hz" 1 / 3 , C 2 = 
5.34 x 10~ 5 Hz -5 / 3 and C 3 = -7.32 x 10~ 9 Hz~ 3 have been obtained by a least-squares fit with 
X 2 = 15.8/12. The v mag profile for GX 340+0 is very similar to those of 4U 1608-52, and Sco X-l 
(see OT99a, Figs. 1-2). 



4. Low Branch Frequency vs kHz QPO Frequency 

The GX 340+0 data of J00 allows one to check the prediction of the model for the low branch 
frequency v^. Taking the profile of Vmagiyvi)^ (Eq. 3), we plot ul m Figure 2. Fixing 5 = 6°. 3, we 
find that our plot for ul and 2vl, calculated according to formula (OT99a, Eq. 9) 

v L = 2u mag (u K / v h ) sin 5, (4) 

fits the data for the observed frequencies of (20-40) Hz and (40-80) Hz respectively. There are three 
observational points in Figure 2 (square points with error bars) which we interpret as the second 
harmonics 2u^. The first harmonics vl for these frequencies are not observed by J00 and they are 
obtained using 2v^ and presented as squares without error bars. The error bars of the observational 
points shown by circle marks are less than the mark size in most cases. It is worth noting that the 
observational points follow the theoretically predicted curve, which shows the signs of saturation 
at high uk- The same type of behavior of vl vs vk is also seen in the Sco X-l data. 



5. Break and Viscous frequencies vs 

Further tests of the TLM can be done using 
break and viscous QPO frequencies vs kHz QPO 



kHz QPO Frequency Correlations 

a comparison of the observed correlations of 
frequencies (J00) with the theoretical depen- 
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dences derived in TLM98 and T099. In Figure 3 (the classification plot), we present the theo- 
retical curves calculated using Eqs. (9, 12-13) in T099 using the dimensionless parameter ax = 
m(a;o/3) 3 / 2 (fo/363 Hz) where vq is the neutron star (NS) spin, m is the NS mass in the solar masses, 
xq is the NS radius in units of Schwarzchild radii. In TO00 (Eq. 4) we approximate the solution 
of Eq (9) of T099 by a polynomial v b = C h P 4 (v K ) for a K = 1.03 and v 0y363 = v /363 Hz = 1. The 
normalization of the theoretical curve is controlled by the constant C& (TO00) which reflects the 
properties of the specific source. For instance, for 4U 1728-34 the observed frequencies are fitted 
by the curve (T099) for which the constant C\ = 1 and Rq = 10 km for 1.4 solar masses (see 
TO00, Eq. 2). It should be pointed out that the NS spin for 4U 1728-34 is known to be vq = 363 
Hz. For GX 340+0 we fit the data points by curves with the parameter vq = 363 Hz for which the 
corresponding constant C& = 7.6. For the source GX 340+0 we found the Lebesgue's measure (for 
a definition, see Titarchuk, Osherovich & Kuznetsov 1999, hereafter TOK) is 16% for ax = 1-03. 
In fact, the theoretical curves are functions of two parameters, ok and vq, and for (a given) fo,363 
that curve with ^o,363 = 1 can be obtained by scaling the argument (see TO00). 

According to solutions presented in T099 and TOK, the break frequency v\> correlates with 
the low Lorentzian frequency vy (which we interpreted as a viscous frequency in the TL model) 
through the power law relation 

Vh oc v v ei . (5) 

Using the data of Wijnands & van der Klis (1999) it was confirmed in TOK that approximately the 
same index is valid for most NS and BH objects. We have also found that the observed correlation 
of v\ } vs vy in GX 340+0 is well described by the same power law dependence (Eq. 5) but with a 
proportionality coefficient 0.08 (which is twice more than that for 4U 1728-34). The dependence of 
vy vs vk is shown in Figure 3 by the red curve. To fit the theoretical curve to the data we use only 
those data points which were specified by J00 as the new discovered frequencies. They pointed out 
that these frequencies well close to half that of the HBO. However looking at the classification in 
Figure 3 we argue that these power spectral features are the same low noise component frequencies - 
in our terms, vy- that were found by Ford and van der Klis (1998) in 4U 1728-34 and by TOK in Sco 
X-l. Furthermore, we have found that at least two power spectra features which were identified by 
J00 as the break frequencies also belong to the viscous branch. The appropriate points are denoted 
by triangle marks in Fig. 3. 

The dependences of z/& and vy on z/k allow us to retrieve information regarding the turbulent 
scale /f p , the sonic velocity v s , the magnetic field strength and structure in the transition layer. We 
discuss these issues in detail in Titarchuk, Bradshaw & Wood (2001). 

Wu (2001) uses the a— prescription - i.e. a ~ {kpVt)/(v 8 H), where H is a half-thickness of 
the disk - to specify the absolute value of the viscous frequency vy. However we argue that the 
introduction of one more free parameter (a) in TL model is not necessary and furthermore it leads 
to very small values of kinematic viscosity, = vtkp/3- One may conclude using Wu's best-fit 
parameters of a{H/R) 2 ~ 10 -3 that vtk p is at least two orders of magnitude less than v r L. On the 
other hand, the Reynold's number 7 is tightly constrained by the observed absolute values of kHz 
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frequencies (TLM98, Fig. 3) and one can find that v t l{ p > Lv r /A. In TLM98 authors have already 
realized this problem with the a— prescription and therefore they have not used it in the model. 

6. Discussion and Conclusions 

In this Letter we have further verified predictions of the transition layer (TL) model for four 
sources by presenting our analysis of Jonker et al.'s GX 340+0 data. Furthermore we demonstrate 
that the right classification of the power spectra features allows us to shed light on their nature 
and origin. The unifying characteristic of TLM is that all QPO features are natural and necessary 
consequence of the adjustment of Keplerian disk to the sub-Keplerian rotation of the neutron star. 
In the transition layer the interaction of the two oscillators may occur because they share the 
common boundary at the outer edge of the transition layer. In this case a strong dependence is 
supposed to exist between all frequencies in the TLM and the fundamental Keplerian frequency z^k 
at the TL outer edge. In our model, assuming high electric conductivity, we describe the frequency 
of the differential rotation of the magnetosphere v mag . We also assume that the NS magnetic field 
has a discontinuity in the equatorial plane similar to the discontinuity of the field in the equatorial 
plane of the solar corona and related current sheet in the heliosphere as described in the model of 
Osherovich, Tzur and Gliner (1984). Both the current sheet and differential rotation are observed 
elements of solar corona and heliosphere (e.g. Hoeksema & Scherrer 1987, Burlaga 1995). 

The TLM model reveals the physical nature of the low branch (HBO) oscillation frequency as 
a low branch of the Keplerian oscillator under the influence of the Coriolis force. The angle 5 as 
a global parameter describes the inclination of the magnetospheric equator to the equatorial plane 
of the disk. Measured locally for different radial distances (therefore different vr)-, <5 may vary 
considerably if the observed oscillations do not correspond to the predicted low Keplerian branch. 
The angle S is calculated using the observed frequencies only. The constancy of 5 shown in Figure 1 
argues in favor of our model. Thus, if a similar rotational splitting effect occurs in BH systems, are 
the needed frequency trio simultaneously observed: a pair of high frequencies of order 200-300 Hz 
(for black hole with mass of 7-10), and a low frequency close to 20 Hz? Strohmayer (2001) reports 
precisely this case for BH candidate GRO J1655-40. He simultaneously observed three frequencies, 
450, 300 and 18 Hz. 

If we assume the rotational splitting paradigm as an origin for these three frequencies, we 
can estimate the rotational frequency of the corona-configuration above the disk Vl/2tt = 167.7 Hz 
using the hybrid relation (OT99, Eq. 7) for = 450 Hz and uk = 300 Hz. Furthermore, if we 
assume that the corona configuration rigidly rotates with the black hole, we can estimate the BH 
dimensionless angular momentum a = cJ/GM 2 = S7GM/c 3 = 3.6 x 10~ 2 . We have also assumed 
a BH mass of about seven solar masses. The low 18 Hz frequency can be interpreted as the low 
branch QPO if 5 = 4.6° (see Eq. 1), which is a typical for NS systems (see Fig. 1). 

Strohmayer (2001) argues that the GRO J1655-40 QPO data, with 0.4 < a < 0.6, could 
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account for the higher frequency, 450 Hz, when identified with the radial epicyclic frequency within 
relativistic disk models (see e.g. Stella, Vietri k Morsink 1999; Psaltis k Norman 2000; Markovich 
2000). He also suggests that the 300 Hz QPO is associated with lower kHz QPOs in neutron 
stars (c./. OT99 and Psaltis et al 1999). However, if this parallel between NS and BH systems 
assumes that the higher frequency QPO twin frequency is due to General relativity (GR) effects, 
the question arises of how this GR model works for NS systems when kHz QPO frequencies of 400- 
500 H (corresponding to 10-12 Schwarzschild radii where GR effects are negligible) are observed. 
Such QPOs, along with higher kHz frequencies, are seen very often in NS systems (see Fig. 4). One 
may conclude that the association between QPOs in NS and BHC along with their interpretation 
within the GR model is somewhat problematic. 

The observed correlations of the viscous and break frequencies versus kHz QPO frequencies 
allow us to retrieve information about the NS spin. We found that for GX 340+0 as well as for Sco 
X-l the NS spin should be close to 360 Hz which is similar to that measured directly in 4U 1728-34 
(Strohmayer et al. 1996). On the other hand there is an indication (TO00) that the NS spin in 4U 
0614+09 can be as high as 700 Hz or more. 

The authors acknowledge discussions with Joe Fainberg, Kent Wood and Reba Bandyopadhyay. 
We are grateful to Peter Jonker for use of his data. We are also grateful to the referee for a thorough 
analysis of the paper content and constructive suggestions 
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Fig. 1. — 5— invariant. Inferred angles 6 between the rotational axis of magnetosphere and the 
normal to the plane of Keplerian oscillations as a function of low peak kHz QPO frequency for four 
sources: GX 340+0 (green), 4U 0614+09 (magenta), Sco X-l (blue), and 4U 1702-42 (red) . The 
5— invariant (eq.[l]) is calculated using the observational frequencies v^, vy. and i>l only. 



Fig. 2.— HBO frequency vs low kHz QPO frequency for GX 340+0 (Jonker et al. 2000). The 
solid line is the theoretical dependence of the low branch oscillation frequencies (black) and the 
second harmonics (blue) on the Keplerian frequency, calculated using multipole expansion of the 
rotational frequency. The points for the first harmonic of the low branch (green squares) are 
obtained by dividing by 2 the observational frequencies of the second harmonic (which are the 
green squares with error bars). 



Fig. 3. — Classification of QPOs in the Z source GX 340+0. Solid lines are the theoretical curves: 
Blue for the Upper hybrid frequency branch, magenta for the second harmonic of the low branch, 
blue for the first harmonic of the low branch, red for the viscous branch and black for the break 
frequency branch. The points with squares for the first harmonic of the low branch are obtained 
by dividing by 2 the observational frequencies of the second harmonic (blue squares) . The observa- 
tional points with black triangles are identified by Jonker et al. (2000) as break frequencies. The 
theoretical curves for viscous and break frequencies branches are constructed using the first six 
viscous frequency points (red asterisks) and the first seven break frequency points (black triangles) 
respectively. 
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